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ABSTRACT 

We present new models for the X-ray photoevaporation of circumstellar discs which suggest that 
the resulting mass loss (occurring mainly over the radial range 10 — 40 AU) may be the dominant 
dispersal mechanism for gas around low mass pre-main sequence stars, contrary to the conclusions of 
previous workers. Our models combine use of the MOCASSIN Monte Carlo radiative transfer code 
and a self-consistent solution of the hydrostatic structure of the irradiated disc. We estimate the 
resulting photoevaporation rates assuming sonic outflow at the surface where the gas temperature 
equals the local escape temperature and derive mass loss rates of ^ 1O~^M0 yr~^, typically a factor 
2 — 10 times lower than the corresponding rates in our previous work (Ercolano et al., 2008) where 
we did not adjust the density structure of the irradiated disc. The somewhat lower rates, and the 
fact that mass loss is concentrated towards slightly smaller radii, result from the puffing up of the 
heated disc at a few AU which partially screens the disc at tens of AU. Our mass loss fluxes agree with 
those of Alexander et al. (2004) but we differ with Alexander et al. in our assessment of the overall 
significance of X-ray photoevaporation, given the large disc radii (and hence emitting area) associated 
with X-ray driven winds. Gorti and Hollenbach (2009), on the other hand, predict considerably lower 
mass loss fluxes than either Alexander et al., (2004) or ourselves and we discuss possible reasons for 
this difference. We highlight the fact that X-ray photoevaporation has two generic advantages for 
disc dispersal compared with photoevaporation by extreme ultraviolet (EUV) photons that are only 
modestly beyond the Lyman limit: the demonstrably large X-ray fluxes of young stars even after 
they have lost their discs and the fact that X-rays are effective at penetrating much larger columns of 
material close to the star. We however stress that our X-ray driven mass loss rates are considerably 
more uncertain than the corresponding rates for EUV photoevaporation (around 10~^'^ Mq yr~^) and 
that this situation will need to be remedied through future radiation hydrodynamical simulations. 
Subject headings: 



1. INTRODUCTION 

The study of the structure and evolution of protoplane- 
tary disks is currently an important area of astronomy as 
it can provide insights into both star and planet forma- 
tion. In recent years there has been considerable debate 
about the final destination of disc material - do discs dis- 
perse through a mixture of accretion onto the star and 
internal consumption by conversion into planets, or are 
they instead dispersed by an extrinsic agent which can 
operate in competition with planet formation? Obvi- 
ously, the answer to this question relates to the lifetime 
of circumstellar discs and the efficiency of planet forma- 
tion around stars of various masses. 

Since disc dispersal by central wind stripping is not 
very efficient (Matsuyama, Johnstone & Hollenbach, 
2009) and since encounters with passing stars is ineffec- 
tive in all but the densest environments (Scally & Clarke 
2001, Pfalzner et al., 2006), and since likewise the effect 
of photoevaporation by neighbouring OB stars requires 
a very rich cluster environment (Johnstone et al., 1998), 
the most likely extrinsic agent of disc dispersal is via 
photoevaporation by energetic (ultraviolet or X-ray) ra- 
diation of the central star. Indeed the detection of Ne ii 



emission in the spectra of e.g Pascucci et al. (2007) and 
Herczeg et al. (2007) is widely interpreted as evidence 
for disc irradiation by high energy photons and poten- 
tially accompanying photoevaporation. Opinion is how- 
ever currently divided as to whether this is best explained 
in terms of irradiation by X-rays (e.g. Glassgold 2007, 
Ercolano 2008b, Meijerink 2008) or ionising ultraviolet 
(EUV) radiation (Gorti & Hollenbach 2008, Alexander 
2008). 

The X-ray properties of low mass pre-main sequence 
stars are well determined (in contrast to the case for the 
EUV radiation from these stars), with flux levels that 
are, if anything, somewhat higher in the case of Weak 
Line (discless) T-Tauri stars than in their disc bearing 
(Classical T Tauri star) counterparts (Preibisch et al., 
2005, Feigelson et al., 2007): although this may result 
from absorption of some of the X-ray emission in accre- 
tion flows close to the star (Gregory et al., 2007, but see 
also Drake et al., 2009), it at least demonstrates that 
the X-rays remain strong throughout and after the disc 
bearing stage and are therefore, in principle, available 
for disc dispersal. The corresponding observational situ- 
ation in the EUV case is much less clear, given the lack of 
suitably high quality ultraviolet spectra, particularly in 
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Weak Line T Tauri stars (Alexander, Clarke & Pringle, 
2005; Kamp & Sammar, 2004). 

Another potential drawback to EUV photoevaporation 
is that, even if one has an accurate estimate of the EUV 
output of the star, it does not necessarily follow that 
these photons can reach the disc, since photons near 
the Lyman limit arc particularly susceptible to absorp- 
tion by even small columns of neutral material close to 
the star. This represents a possible counter- argument to 
Ne II emission being excited by EUV radiation in most 
cases. It is usually therefore argued (e.g. Gorti & Hol- 
lenbach, 2009) that EUV photoevaporation cannot be 
effective until late times when the subsiding of both ac- 
cretion and outflow activity should allow the disc to be 
exposed to ionising photons. From this point onwards, 
the theory is well developed (Hollenbach et al., 1994, 
Clarke et al., 2001, Alexander et al., 2006a, b), with a 
predicted photoevaporation rate of 10^^°Mq yr^^, for 
an EUV luminosity of ~10^^s~^ (the rate goes as the 
square root of the assumed luminosity), and the inter- 
play between photoevaporation and accretion producing 
a characteristic evolutionary sequence involving the cre- 
ation and rapid outward expansion of an inner disc hole. 
Thus EUV photoevaporation provides an alternative to 
planet formation as the origin of inner hole (transition 
disc) systems (see e.g. Najita et al., 2007). 

Notwithstanding this success of EUV photoevapora- 
tion at late times in removing the last few Jupiter masses 
of gas from the disc, there has also been some discus- 
sion of whether other radiation sources, such as X-rays 
or non-ionising (FUV) ultraviolet photons, might not 
be much more effective disc photoevaporation agents at 
early times and might indeed be capable of removing 
much larger quantities of disc gas. Gorti & Hollenbach 
(2009) investigated this issue and argued that FUV pho- 
toevaporation is indeed a very significant effect at early 
times when accretion produces a strong FUV flux: their 
predicted FUV mass loss rates of > 1O~^M0 yr~^ exceed 
EUV rates by more than two orders of magnitude. They 
however concluded, by contrast, that X-rays, in the ab- 
sence of FUV radiation, are considerably less important 
even than the EUV, a conclusion also shared by Alexan- 
der et al. (2004). 

Another recent study of X-ray photoevaporation was 
provided by Ercolano et al 2008b, henceforth Paper I. 
While Paper I mainly aimed at understanding the ther- 
mochemical structure of the irradiated regions of the disk 
and at identifying gas-phase emission line diagnostics 
(see also Ercolano, Drake & Clarke 2008c), it also yielded 
an estimate for the X-ray photoevaporation rate that was 
considerably higher than those found by Gorti & Hollen- 
bach (2009) and Alexander et al (2004), being around 
1O"*^M0 yr^i. This estimate was however very prelimi- 
nary, due not only to the general difficulty of estimating 
mass loss rates from static models (see Section 6), but on 
account of the fact that the disc's density structure was 
not modified in response to heating of its upper layers. 
Instead it was assumed that the disc adopted the fixed 
structure that it would have in the case that its tempera- 
ture structure was controlled by full thermal coupling be- 
tween dust and gas (d'Alessio et al., 1998). This is clearly 
inconsistent given that the gas and dust temperatures in 
the X-ray irradiated disc turn out to be completely de- 
coupled down to a column depth of ~5xl0^^cm~^. In 



contrast. Gorti & Hollenbach (2009) and Alexander et al. 
(2004) both iterated to a self-consistent hydrostatic den- 
sity structure for the irradiated disc. Since these authors 
had previously concluded that X-ray photoevaporation is 
a minor effect, it would appear necessary to check how a 
similar iteration onto a hydrostatic structure would affect 
the initial estimate of Paper I. The reason for pursuing 
a study that is broadly similar to those of Gorti & Hol- 
lenbach (2009) and Alexander et al. (2004), is that there 
are a number of differences in the way that the radiative 
transfer is treated in the three works. We employ a fully 
three dimensional Monte Carlo approach to the radia- 
tive transfer through dust and gas, and, while - in the 
present paper - we iterate on the density profile at each 
cylindrical radius, the final radiation field and temper- 
ature structure represents the thermal equilibrium con- 
dition for a given density field, without any simplifying 
assumptions about the geometry of the resulting radia- 
tive fiux vectors. The other two works instead employ 
methods that are pseudo-one dimensional in that they 
integrate along a range of lines of sight from the X-ray 
soiirce. There are also some other differences concerning 
the assumed spectrum of the X-ray source and some de- 
tails about the detailed interaction between X-rays and 
gas that are discussed in more detail in Section 5. 

The calculations presented here thus aim specifically at 
providing an improved estimate of X-ray photoevapora- 
tion rates, together with a comparison with the results of 
using irradiating spectra that also include an EUV com- 
ponent. For computational tractability, only the inner 
50AU of the disk, from which significant photoevapora- 
tion is found, are considered here and therefore a full 
emission line spectrum is not available. This will be in- 
vestigated in future work. 

A description of the computational methods is given 
in Section 2. Our improved photoevaporation rates are 
presented and discussed in Section 3, while Section 4 
describes the resulting thermal and ionisation structures. 
Section 5 analyses differences in methods and results 
compared with previous studies and Section 6 sets out 
the main uncertainties involved in estimating mass loss 
rates from hydrostatic models. Section 7 summarises our 
main results. 

2. PHOTOIONISATION MODELS IN 
HYDROSTATIC EQUILIBRIUM 

We have used a modified version of the 3D photoion- 
isation and dust radiative transfer code MOCASSIN (Er- 
colano et al., 2003; 2005; 2008a). This code uses a Monte 
Carlo approach to the three-dimensional transfer of radi- 
ation, allowing for the self-consistent transfer of both the 
stellar (primary) and diffuse (secondary) components of 
the radiation field. Details of atomic data and physical 
processes used for the photoionisation calculations are 
given by Ercolano et al. (2003, 2008a) and Ercolano & 
Storey (2006). The main changes to the thermal balance 
and dimensional setup of the code have been described 
in detail in Paper I and are briefly summarised here, 
(i) A treatment of viscous heating is included using the 
formulation for a thin disk by Pringle (1981). (ii) Gas 
cooling by collisions of grains with a mixture of atomic 
and molecular hydrogen is calculated using the method 
of Hollenbach & McKee (1979), also discussed by Glass- 
gold et al. (2004). (iii) mocassin was adapted to allow 
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for two dimensional simulations in order to exploit the 
symmetry of the system. 

In Paper I, viscous accretion heating was included us- 
ing the standard treatment for a thin disk (Pringle 1981) 
with viscosity parameter a =0.01. This corresponds to 
an accretion rate of the order of 10'^^ MQ/yr. This is 
rather high, in particular as we are interested in the ef- 
fects of X-ray photoevaporation in later stages of disk 
evolution when accretion has turned off. This is achieved 
in our models by setting a =10"^" which effectively ren- 
ders accretion heating negligible. We have also improved 
the model in Paper I by including excitation of C i and 
O I fine structure lines by collisions with neutral hydro- 
gen (only collisions with electrons were previously in- 
cluded, see Ercolano, Drake & Clarke, 2008 for more de- 
tails) . 

2.1. Hydrostatic Equilibrium 

We have further modified the MOCASSIN code to in- 
clude a self-consistent calculation of the ID vertical 
hydrostatic equilibrium structure (HSE) of the X-ray 
heated disk. We used the methods described by Alexan- 
der, Clarke and Pringle (2004, ACP04) to update the 
disk density structure after convergence of the thermal 
and ionisation structure based on the previous HSE iter- 
ation. 

Our initial density distribution is the same as that used 
in Paper I and it consists of the disk model structure 
calculated by D'Alessio et al. (1998), selected to be that 
which best fits the median SED of T Tauri stars in Taurus 
(d'Alessio, 2003). This model was calculated under the 
assumption of full thermal coupling between dust and gas 
and featured an irradiating star of mass 0.7 Mq, radius 
2.5 Rq, and an effective temperature of 4000 K. Disk 
parameters include a mass accretion rate of 10~^ Mq 
yr~^ and a viscosity parameter a = 0.01. The total mass 
of the disk is 0.027 M© and the outer radius of --500 AU. 
The surface density in the disk is inversely proportional 
to the cylindrical radius. 

We perform HSE calculations in the vertical direction, 
using as a lower boundary the point on the z-axis where 
the gas and dust temperatures become decoupled. As in 
Paper I, we keep the dust temperatures fixed to the val- 
ues calculated by D'Alessio et al. (2001) given that with- 
out treating the stellar or interstellar optical and far-UV 
fields, the dust temperatures we would obtain from our 
model would be incorrect. Calculating the hydrostatic 
equilibrium from the decoupling point has the advantage 
of keeping a realistic density structure below this point 
even when only X-ray irradiation is considered (i.e. not 
treating the photospheric emission of the pre-main se- 
quence star). 

The density structure from the lower boundary point 
is calculated using a slightly modified version of equa- 
tion A9 of ACP04, as shown below. At a given cylindri- 
cal radius, the number density of the g clS 3tS 3j function 
of height above the midplane, z, is given by 



atom and k is the Boltzmann constant. (3{z) is calculated 
as follows 
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where and Tfj are density and temperature at the 
lower boundary, is the mass of the star, G is the 
gravitational constant, mn is the mass of the hydrogen 



fi(z')z' 



^ T{z'){z'^ + R^)t 



-dz' 



(2) 



where ii{z) is the mean molecular weight of the gas at z 
and R is the cylindrical radius. 

Equations [1] and [2] differ from equation A9 of ACP04 
in that we do not assume z << R. We note, however, 
that a full treatment would consist of a two-dimensional 
HSE in the radial and vertical directions - with gravity 
balanced by pressure in the vertical direction and with 
gravity balanced by a combination of pressure and cen- 
trifugal forces in the radial directions. In the case of an 
irradiated disk with significant vertical extent the pres- 
sure terms are very complicated and a full solution is not 
attempted in this paper. It is indeed unclear whether 
a fully two-dimensional hydrostatic solution would im- 
prove our calculations, particularly in the upper layers 
of the disk, where the assumption of hydrostatic equilib- 
rium is altogether poor given that material is likely to 
be flowing in a photoevaporative wind. However a full 
hydrodynamical flow solution is beyond the scope of this 
paper. 

2.2. The dust model 

Although we do not calculate the dust temperatures, 
the effects of gas-dust interactions on the temperature 
structure of the gas, are still taken into account in our gas 
thermal balance as well as in the radiative transfer, where 
the competition between dust and gas for the absorp- 
tion of radiation is properly treated (see Ercolano et al. , 
2005). Our dust treatment consists of a standard MRN- 
type model (Mathis, Rumpl & Nordsiek, 1977) with 
minimum and maximum grain radii Omin — 0.005/im 
and Omax = 0.25/x. The disk structure of d'Alessio et 
al. (1998), which provides our starting density distri- 
bution, employs a bimodal dust distribution, where at- 
mospheric dust follows the standard MRN model and 
interior dust consists of larger grains with a size distri- 
bution still described by a power law of index of —3.5, 
but with Omin — 0.005/im and Omax = 1mm. The transi- 
tion between atmospheric and interior dust occurs at an 
height of 0.1 times the midplane gas scale height. As we 
are mainly interested in regions well above the transition 
point, we have chosen for simplicity to use atmospheric 
dust everywhere. Following d'Alessio et al. (2001) the 
dust to gas mass ratio of graphite is 0.00025 and that of 
siHcates is 0.0004. 

High energy dust absorption and scattering coefficients 
are calculated from the dielectric constants for graphite 
and silicates of Laor & Draine (1993), which extend 
to the X-ray domain. Spherical grains are assumed 
and we use the standard Mie scattering series expan- 
sion for complex refractive indices, x\m\ < 1000, where 
X = 2 • a/ A is the scattering parameter (see Laor & 
Draine, 1993). For x\m\ > 1000 and x\m - 1| < 0.001 
we use Rayleigh-Gans theory (Bohren & Huffman, 1983), 
and for x\m\ > 1000 and x\m — 1| > 0.001 we use the 
treatment specified by Laor & Draine (1993), which is 
based on geometric optics. 

2.3. T- Tauri EUV- X-ray Spectra 
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Fig. 1. — Isothermal model spectrum for Log(Tx) = 7.2. 

We illuminate the disk using synthetic X-ray spectra 
representative of the high-energy emission from a typical 
T Tauri star corona. This coronal X-ray emission arises 
from a collision-dominated optically-thin plasma. 

We first assumed a single X-ray temperature, 
\og(Tx) — 7.2, and a number of luminosities spanning 
2 orders of magnitudes around our base luminosity of 
Lx(0.1-10 keV)== 2 x lO^" erg/sec (Model XSOHlOLxl). 
The labelling scheme of our models is explained in Sec- 
tion 2.4. The bulk of the luminosity in these models 
is emitted in the 1-2 keV region as shown in Figure [TJ 
The same unabsorbed isothermal model spectrum for 
Log(Tx) = 7.2 was also used in Paper I and we refer 
to that work for a more detailed description of the com- 
putational methods and atomic data used. 

While such a model provides an accurate description 
of the X-ray spectrum, it represents only a lower limit to 
the unabsorbed flux at EUV wavelengths 100-912 A) 
where cooler plasma can contribute significant flux. In 
Paper I we argued that for most of a disk's lifetime EUV 
photons would not easily reach the disk due to absorp- 
tion from circumstellar material and restricted our in- 
vestigation to the X-ray dominated irradiation spectrum 
above. Here we expand our study to explore the effects 
that different amounts of EUV irradiation would have 
on a disk that is also simultaneously irradiated by X- 
ray. In a companion paper (Glassgold et al., 2009 in 
prep.) we study the relative importance of EUV and X- 
rays on ionisation and heating rates, while here we focus 
on their effects on photoevaporation. To this aim, we 
employ a more sophisticated multi-temperature plasma 
model developed to provide an accurate representation 
of the ionizing coronal spectrum throughout the EUV- 
X-ray range. We have not included any contributions to 
the ionizing flux from accretion here; this will be consid- 
ered in future work. 

Observed line-of-sight absorbing columns toward 
T Tauri stars are typically greater than 10^° atoms cm~^ 
(e.g. Feigelson & Montmerle 1999), and consequently 
their EUV spectra are completely obscured from our view 
by photoelectric absorption. While this would appear to 
be a rather serious impediment to producing a realistic 
T Tauri EUV spectrum, there are much more nearby 
and relatively unabsorbed sources that provide a tem- 
plate: the RS CVn class of active binaries that comprise 
at least one evolved star lying on the sub-giant or gi- 
ant branch. The coronal emission of RS CVn binaries 
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Fig. 2. — Emission measure distributions as a function of tem- 
perature. 

is largely dominated by the evolved component that is 
likely to be magnetically more similar to a T Tauri star 
than are active main-sequence stars with shallower con- 
vection zones. Rotation periods of RS CVn binaries are 
also similar to those of T Tauri stars and typically lie in 
the range 1-10 days (e.g., Strassmeier et al. 1993, Get- 
man et al., 2005). 

The expected resemblance of the coronal properties of 
RS CVn-type binaries and T Tauri stars is borne out 
in X-ray studies. For example, the Einstein survey of 
late-type stars by Schmitt et al. (1990) obtained plasma 
temperatures for RS CVns in the range 7 < logT < 7.8 
with a mean of order logT 7.3, and X-ray luminosities 
in the range lO'^'^ < Lx ^ lO'^^'^ erg s~^. These ranges 
are very similar to those found for T Tauri stars of the 
Orion Nebular Cluster (ONC) in the mass range 1-2Mq 
by Preibisch et al. (2005). On this basis, the use of 
RS CVn EUV spectra as a proxy for the coronal spectra 
of T Tauri stars seems well-justified. 

Detailed studies of the outer atmospheric temperature 
structure of RS CVn-type binaries also finds similarity in 
the transition regions and lower coronac of different stars. 
Sanz-Forcada et al. (2002) studied four active RS CVns 
in detail from UV to the shortest EUV wavelengths using 
spectra obtained from the International Ultraviolet Ex- 
plorer (lUE), the Extreme Ultraviolet Explorer (EUVE) 
and the Orbiting and Retrievable Far and Extreme Ultra- 
violet Spectrometer (ORFEUS). They derived emission 
measure distributions as a function of temperature for 
the range logT = 4.5-7.2 and we use a representation of 
these for the emission measure distribution adopted as 
the basis of our T Tauri model. This emission measure 
distribution is illustrated in Figure [2l The Sanz-Forcada 
et al. (2002) emission measures are rather ill-constrained 
for temperatures log T > 7 owing to the paucity of spec- 
tral line diagnostics in the EUV range at these tempera- 
tures. For this range we have adopted a smooth function 
that peaks at logT = 7.25 and drops off fairly sharply for 
temperatures logT > 7.5. The temperature of the peak 
is based on the two-temperature optically-thin thermal 
plasma model fits to the large sample of Chandra ONC 
T Tauri spectra by Maggio et al. (2007). 

Coronal spectra based on the emission measure distri- 
bution in Figure[2]were computed as described in Paper I 
using the PINTofALE IDL software suite-^ (Kashyap 

1 PINTofALE is freely available from 
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Fig. 3. — Model coronal spectra based on the emission measure 
distribution shown in Figure[2] The black line represents the unat- 
tenuated spectrum, while coloured lines refer to spectra transmit- 
ted through semi-neutral columns of the values shown. For clarity, 
the spectra are shown in arbitrary units and are not normalised to 
the X-ray luminosity values used in the modelling. 

& Drake 2000), using line and continuum emissivities 
from the CHIANTI compilation of atomic data (Landi 
& Phillips 2006, and references therein), ion populations 
from Mazzotta et al. (1998), and the solar chemical com- 
position of Grevesse & Sauval (1998). Figure [3] shows the 
resulting EUV-I- X-ray spectrum produced (black Hue); 
the emergent luminosity is roughly the same in the EUV 
(13.6< V <100 eV) and X-ray regions {v > lOOeV). 

In order to account for absorption of coronal EUV 
and X-ray photons in the immediate local stellar and 
circumstellar environment, disk irradiation calculations 
were also performed for spectra attenuated by photoelec- 
tric absorption corresponding to neutral hydrogen col- 
umn densities in the range Nh = 10^^-10^^ cm^. The 
gas transmittance was computed using the cross-sections 
compiled by Wilms, Allen, & McCray (2000) as imple- 
mented in PINTofALE. The effects of different 'pre- 
screening' columns on the radiation impinging on the 
disk are shown by the coloured lines in Figure [31 where 
the black, red, green, blue, magenta and cyan lines show 
the transmitted spectrum through columns of 0., 10^^, 
10^^, 10^°, 10^^ and 10^^ cm~^, respectively. The corre- 
sponding models are labeled FS0H2, FS18H2, FS19H2, 
FS20H2, FS21H2, FS22H2. For the disk photoionisation 
calculations, spectra were normalised to an X-ray lumi- 
nosity Lx(0.1 - 10 keV) = 2 • 10^°erg/sec. 

Finally, in Paper I, we assumed the height of the stellar 
irradiating source to be 10 stellar radii (R*) above the 
disk mid-plane, as was adopted in earlier studies (e.g. 
Glassgold et al, 1997, Igea & Glassgold 1999) This repre- 
sents quite a large scale height for coronal emission, and 
the source of quiescent X-rays is likely to reside much 
closer to the star. Indeed, Flaccomio et al. (2005) found 
evidence for rotationally modulated X-rays from T Tauri 
stars in the Orion Nebular Cluster, indicating that the 
dominant emission from those stars lies at scale heights of 
order the stellar radius or below. Gregory et al. (2006) 
also find that simple extrapolations of T Tauri surface 
magnetograms suggest compact coronal structures. We 
therefore adopt here a height of 2 R* for the X-ray source, 
corresponding to the X-ray emission originating from a 
point located 1 R* above the stellar pole. We also in- 

Ihttp: / /hea-www. harvard.edu/PINTofALE] 



vestigate the effect of this scale height choice and find 
that changing the source height between 2 and 10 R* 
has relatively little influence on the results. 

2.4. Model specification 

Table [1] contains a summary of the parameters used 
for the models presented in this paper. In general mod- 
els are labeled as nSsHhLxl, where 'n' is 'F' for models 
using EUV-f X-ray irradiating spectra (Figure [S]) and 'X' 
for models using the single temperature Log(Tx) — 7.2 
spectrum, 's' indicates the pre-screening column and it is 
for no pre-screening and 22 for a pre-screening column 
of lO^^cm"^. 'h' indicates the height of the irradiating 
source which is 2 or 10 (stellar radii). '1' indicates the 
X-ray luminosity used, e.g. 01, 1, 10 for Lx = 2 ■ 10^^, 
2-l0^°, 2-103\ respectively. 

We have adopted the following elemental abundances, 
given as number densities with respect to hydrogen, 
He/H = 0.1, C/H = 1.4x10-", N/H = 8.32xl0-^ 
O/H = 3.2xl0-*, Ne/H= 1.2x10-4, Mg/H = 1.1 x 10"^ 
Si/H = 1.7x10-^ S/H = 2.8x10-^ These are solar 
abundances (Asplund et al. 2005) depleted according to 
Savage & Sembach (1996). 

Given the differences in the code and input parameters 
from those used in Paper I, we also produced a set of 
fixed density structure models in order to compare with 
the HSE and isolate the effects of the latter. This is made 
necessary particularly by the fact that the HSE models 
are likely to overestimate the densities in the upper lay- 
ers (which would be in reality part of a photoevaporative 
flow), hence fixed density structure models, which under- 
estimate the density of the upper layers, can be used to 
provide an educated guess of the uncertainties in the re- 
sults. "Fixed grid" models are listed in Tabled] with the 
suffix 'FG'. 

3. X-RAY PHOTOEVAPORATION RATES AND 
TIMESCALES FOR DISK DISPERSAL 

We estimate photoevaporation rates assuming that the 
mass loss rate per unit surface is roughly S = pCs, where 
p and Cs are, respectively, the gas density and the sound 
speed evaluated at the base of the flow. The base of the 
flow is defined as the location in the disk atmosphere 
where the gas temperature exceeds the local escape tem- 
perature, Tes = GmHM^,/kR, where M^, is the stellar 
mass and R is the radial distance. The total mass loss 
rates M are obtained by integrating E over the disk sur- 
face. This simple approach, as discussed in Paper I, 
introduces significant uncertainties in the calculation of 
mass loss rates. These are further discussed in Section [5] 
and[S] 

The last column of Tabic [T] lists the mass loss rates 
obtained by our models. The sensitivity of these rates to 
the model parameters are discussed in more detail below. 
Our general results show that soft X-rays (roughly a few 
100 eV) are most efficient at driving a photoevaporative 
wind, producing a mass loss rate of the order of 10~'^ 
Me/yr. 

Illuminating spectrum and pre-screening columns — The 
values in Table [1] show that both for HSE and fixed den- 
sity structure models the total mass loss rates are not 
very sensitive to the inclusion of EUV (<100eV) radi- 
ation. The models show very similar values of M for 
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TABLE 1 

Model input parameters and predicted mass loss rates from photoevaporation. 'Spectrum' refers to the irradiating 

SPECTRUM, see TEXT IN SECTION 2.3. 'COLUMN' IS THE COLUMN DENSITY OF CIRCUMSTELLAR SCREENING MATERIAL (PRE-SCREENING 

column). 'Height' is the height of the irradiating source. Lx is the X-ray luminosity defined from 100 eV < E < 10 keV. 
'Hydro' refers to whether the models are in hydrostatic equilibrium (TRUE) or a fixed density structure was used 

(FALSE) 





Spectrum 


Column 


Height 


Lx 


Hydro 


M 








[cm~^] 


[R*] 


[erg /sec] 




[Mo/yr] 


XS0H2Lxl 




X-ray 


0. 


2 


2.e30 


TRUE 


2.0e-9 


XSOHlOLxl 




X-ray 


0. 


10 


2.e30 


TRUE 


1.2e-9 


XSOHlOLxOl 




X-ray 


0. 


10 


2.e29 


TRUE 


2.4e-10 


XS0H10Lx02 




X-ray 


0. 


10 


4.e29 


TRUE 


4.5e-10 


XS0H10Lx04 




X-ray 


0. 


10 


8.e29 


TRUE 


7.0e-10 


XSOHlOLxOS 




X-ray 


0. 


10 


1.6e30 


TRUE 


l.le-9 


XS0H10Lx2 




X-ray 


0. 


10 


4.e30 


TRUE 


2.2e-9 


XS0H10Lx4 




X-ray 


0. 


10 


8.e31 


TRUE 


4.0e-9 


XSOHlOLxlO 




X-ray 


0. 


10 


2.e31 


TRUE 


5.9e-9 


XS0H10Lx20 




X-ray 


0. 


10 


4.e31 


TRUE 


l.le-8 


FSOHlOLxl 


X- 


ray-(-EUV 


0. 


10 


2.e30 


TRUE 


3.5e-9 


FS0H2Lxl 


X- 


ray+EUV 


0. 


2 


2.e30 


TRUE 


4.5e-9 


FS18H2Lxl 


X- 


ray-(-EUV 


1018 


2 


2.e30 


TRUE 


4.5e-9 


FS19H2Lxl 


X- 


ray-(-EUV 


10« 


2 


2.e30 


TRUE 


4.2e-9 


FS20H2Lxl 


X- 


ray-l-EUV 


1020 


2 


2.e30 


TRUE 


4.0e-9 


FS21H2Lxl 


X- 


ray+EUV 


1021 


2 


2.e30 


TRUE 


2.7e-ll 


FS22H2Lxl 


X- 


ray+EUV 


1022 


2 


2.e30 


TRUE 




XS0H2Lxl.FG 




X-ray 


0. 


2 


2.C30 


FALSE 


l.le-8 


XSOHlOLxLFG 




X-ray 


0. 


10 


2.C30 


FALSE 


8.1e-9 


FSOHlOLxLFG 


X- 


ray+EUV 


0. 


10 


2.e30 


FALSE 


2.2e-8 


FS0H2LX.FG 


X- 


ray+EUV 


0. 


2 


2.e30 


FALSE 


2.2e-8 


FS18H2Lxl.FG 


X- 


ray+EUV 


1018 


2 


2.e30 


FALSE 


1.5e-8 


FS19H2Lxl.FG 


X- 


ray+EUV 


10" 


2 


2.e30 


FALSE 


1.3e-8 


FS20H2Lxl.FG 


X- 


ray+EUV 


1020 


2 


2.e30 


FALSE 


7.7e-9 


FS21H2Lxl_FG 


X- 


ray+EUV 


1021 


2 


2.e30 


FALSE 


1.3e-9 


FS22H2Lxl.FG 


X- 


ray+EUV 


1022 


2 


2.C30 


FALSE 






□ 10 20 30 40 50 10 20 30 40 50 

Rodius [AU] Rodius [AU] 



Fig. 4. — Surface mass loss rates (left) and total mass loss rates 
(right) of hydrostatic equilibrium X-ray+EUV models with differ- 
ent pre-screening columns as indicated by the model labels. The 
solid lines correspond to models irradiated by a source located at a 
height of 2 R* , while the dashed line corresponds to an unscreened 
X-ray+EUV model irradiated by a source at 10 R, . 

pre-screening columns ranging from 0. to 10^*^0™^^, im- 
plying that the radiation screened by columns of up to 
10^°CTO~^ (energies less than lOOeV) is less efheient than 
higher energy radiation at dispersing the disk by photo- 
evaporation. This can be understood by examining the 
surface mass loss rates, S, and total mass loss rates, M, 
of models with different pre-screening columns, as shown 
for the hydrostatic equilibrium models in Figure [H Al- 
though E at 5-lOAU decreases as the pre-screening col- 
umn increases (see left panel of Figure [4]) , however E at 



larger radii (20-30 AU) is roughly invariant for columns 
up to lO^^cm"^. Total mass loss rates are dominated 
by the surface emission at larger radii, due to the ra- 
dius squared dependence of the mass loss rate integral, 
and therefore they also appear to be invariant for absorb- 
ing columns up to 10^'^ cm~^. This is shown in the right 
panel of Figured] where the quantity 2ttR^'E, is plotted as 
a function of radius, R. We plot 2ttR'^'S as this quantity 
is generally used as a rough estimate of M; in Table [1] 
however, we give the values of M obtained by formally 
integrating E over the whole disk. 

The reason why S at radii in the 20-30 AU range ap- 
pears to be roughly the same for models with obscuring 
columns up to lO^^cm"^ (which screen out EUV radia- 
tion) is that the thermal structure of the disk at this ra- 
dial distance is hardly affected by EUV radiation which 
is very weak at these radii due to large neutral H opacity 
of the disk at EUV wavelengths. Radiation at soft X-ray 
wavelengths, on the other hand, is optimal for heating 
this region of the disk. 

For columns larger than lO^'-'cm'^^ we see a dras- 
tic drop in total mass loss rates, as shown in Ta- 
ble 1. A model with a pre-screening column of lO^^cm"^ 
(FS21H2Lxl; cyan line in Figure |4]) is approximately 
a factor of 200 lower than for models with lower pre- 
screening columns and a model with a lO^^cm"^ pre- 
screenmg column (FS22H2Lxl; not plotted in Figure H]) 
shows no photoevaporation at all. This is due to the 
fact that these high pre-screening column models, while 
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Rodlus [AU] Rodius [AU] 

Fig. 5. — Surface mass loss rates (left) and total mass loss 
rates (right) of fixed density X-ray+EUV models with different pre- 
screening columns as indicated by the model labels. The solid linos 
correspond to models irradiated by a source located at a height 
of 2 R*, while the dashed line corresponds to an unscreened X- 
ray+EUV model irradiated by a source at 10 R» . 
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Fig. 6. — Surface mass loss rates (left) and total mass loss rates 
(right) of hydrostatic equilibrium single-temperature X-ray mod- 
els of different luminosities as indicated by the model labels. All 
models were irradiated by a source located at a height of 10 R» . 

having the same total X-ray luminosity as less screened 
ones, have a much harder spectrum (see Figure [3l but 
note that for clarity the spectra shown in this figure have 
not been normalised). The reason for keeping the total 
X-ray luminosity constant is to be able to isolate the ef- 
fects of the irradiating spectral shape. "Hard" X-rays (> 
1 keV) have very large penetration depths and therefore 
they do not heat the upper layers of the disk as efficiently 
as their lower energy counterparts. Their energy is "dis- 
persed" in a larger column and, while these photons are 
vital to achieve a low-level ionisation deep in the disk 
interior, which is a basic requirement of the magneto ro- 
tational instability theory (Balbus and Hawley, 1991), 
their contribution to the heating of the upper disk layers 
and therefore to disk dispersal is insignificant. We note 
at this point that previous work by Gorti & HoUenbach 
(2009, GH09), who used a rather hard X-ray spectrum 
to irradiate their disks also obtained very low photoevap- 
oration rates (see discussion in Section [5]). 

Height of the illuminating source — Results of models with 
the illuminating source at a height H = 2 i?* are not 
too dissimilar from those with H = 10 i?* (as assumed 
in Paper I). The radial dependences of S and M for a 
model with H = 10 i?* (FSOHlOLxl) are represented by 



the dashed lines in the left and right panels of Figure |4l 
respectively. S at 5-10 AU for the H = 10 i?, model 
(FSOHlOLxl) is comparable to the values obtained for 
the H = 2 model (FS0H2Lxl), at -20 AU, however 
the H = 10 i?* model produces about a factor of two 
higher S, but negligible levels of photoevaporation at 
larger radii (30-40 AU). This results in E x 2ttR'^ peak- 
ing at —20 AU rather than ~35 AU as for the model 
with H = 2 i?* . The reason for this inward shift of the 
peak mass loss flow is due to the fact that radiation emit- 
ted from H = 10 i?* hits the inner disk traveling along 
rays that are closer to the normal to the disk surface, 
resulting in a more efficient heating of these regions. A 
warmer inner disk is more puffed up and more likely to 
screen the material at larger radii, hence producing lower 
photoevaporation rates than in the H = 2 i?* case. 

The above interpretation is borne out by the fixed den- 
sity structure model results for H = 10 i?* and H = 2 i?* 
(FS0H10Lxl_FG and FS0H2Lxl.FG). Here the density 
distribution of the disk is fixed and equal for both models 
and, indeed, we see no difference in the mass loss rates 
(see Table [T] and dashed lines in Figure [5]). 

We finally note that the total mass loss rate obtained 
from the HSE model with H 10 i?, (FSOHlOLxl) is 
only about 20% smaller than that obtained for the same 
model with with H = 2 i?* FSOHlOLxl (FS0H2Lxl), and 
therefore not a real cause for concern, given the (larger) 
error implicit to our pCg estimate of E (see Section 6). 

X-ray luminosity — Since X-ray luminosities of T Tauri 
stars range over two orders of magnitude or more, it is 
of interest to investigate the general trends of photoe- 
vaporative mass loss rates with varying model X-ray lu- 
minosity. For a discussion of the connection between 
photoevaporation, accretion rate and X-ray luminosities 
we refer to a companion paper by Drake et al. (2009). 

To eliminate the effects of EUV radiation, which may 
complicate any trends due to X-rays, we irradiate the 
disk with the single temperature X-ray model already 
used in Paper I (Figure [T]) . Figure [6] shows the results 
for the radial distribution of E and M for models with 
X-ray luminosities varying between 2 • lO^^erg/ sec and 
4 • lO^^erg/sec. The total mass loss rates are plotted 
as a function of X-ray luminosity in Figure llOi which 
shows a roughly linear correlation of these quantities. 
The increase in X-ray luminosity means that deeper lay- 
ers (higher density) in the disk can be heated to tempera- 
tures in excess of the local escape temperatures, directly 
resulting in an increase in the mass loss rates. In addition 
the higher temperatures cause a vertical expansion of the 
heated surface material, hence increasing the solid angle 
directly irradiated by the source. However, the puffed up 
inner disk also obscures the disk at higher radii, hence 
inhibiting further increases in mass loss rate. This ex- 
plains the slightly weaker than linear Lx-M relation in 
Figure [TUl for high Lx- 

Dtsk structure and Hydrostatic Equilibrium — The results 
of models that use the HSE procedure described in Sec- 
tion 12.11 are compared here with those obtained using a 
fixed density structure for the same parameters. 

A casual inspection of Table 1 shows that hydrostatic 
equilibrium models produce smaller mass loss rates than 
the fixed density structure models with the same param- 
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Fig. 7. — Gas density distribution for the unscreened X- 
ray-l-EUV model in HSE (FSOHOLxl). The thick black line in- 
dicates the location of the photoevaporative flow base. 




-1.0 -0.5 0.0 0.5 1.0 1.5 

log Radius [AU] 



Fig. 8. — Gas density distribution for the unscreened X- 
ray-l-EUV model with fixed grids (FSOHOLxl.FG). The thick black 
line indicates the location of the photoevaporative flow base. 

eters. As in the case of the height of the illuminating 
source, the difference is due to the fact that, unlike fixed 
density structure models, heated material in the disk of 
HSE models puffs up and screens material at larger radii 
from the source, producing lower photoevaporation rates. 

Figures [7] and [H] show the two-dimensional density dis- 
tribution of the disk for unscreened X-ray+EUV models 
respectively in HSE and with a fixed density (FS0H2Lxl 
and FS0H2Lxl_FG). The base of the photoevaporative 
envelope is marked by the thick black line. The physical 
height, density and temperatures of the flow surfaces are 
shown in Figure [9l 

Figure shows the radial distribution of S and M for 
the fixed density structure models. A comparison with 
Figured! which shows the same quantities for the HSE 
models, promptly reveals that, unlike fixed density struc- 
ture models, HSE models show photoevaporation from 
radii smaller than ^^10 AU. This is due to the fact that 
heated gas in HSE models is allowed to expand in the 
z-direction, moving away from the star to locations of 
lower gravitational potential where escape temperatures 
are lower, thus establishing a photoevaporative flow. 

Furthermore gas at larger radii in fixed density struc- 
ture models does not suffer from attenuation from puffed 
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Fig. 9. — Height, gas density and temperature of the photoe- 
vaporative flow surface for the unscreened X-ray-|-EUV models in 
HSE (black asterisks, FS0H2Lxl) and with fixed density (red tri- 
angles, FS0H2Lx1_FG). At small cylindrical radii (R<10AU) the 
temperatures required for the gas to escape are only reached at 
large heights above the midplane in the hydrostatic equilibrium 
models. The fixed grid models show no photoevaporation at all at 
R<10AU due to the fact that the disk atmosphere is not allowed 
to expand and the escape temperatures are therefore never reached 
by the gas in these regions (see Section 3). 

up inner material, resulting in substantial mass loss rates 
being supported in a larger portion of the disk. For ease 
of computation, our simulations only included the inner 
50AU of the disk. As it is apparent from the left panel of 
Figure m this is not sufficient to estimate total mass loss 
rates from low prc-scrcening (< 10^^ cm^^) fixed density 
structure models, and the values of M given in Table 1 
should be regarded as lower limits. 

4. DENSITY AND TEMPERATURE STRUCTURE 

A detailed discussion of the thermochemical and den- 
sity structures predicted by our models will be presented 
in future work, where predictions for gas emission line 
diagnostics will also be compiled for a larger disk model. 
Nonetheless, it is helpful to illustrate briefly here the ef- 
fects of X-ray and EUV irradiation on the density and 
temperature distribution of gas in the disk, in particu- 
lar for those regions where the photoevaporative flow is 
predicted to be significant. 

Figure [Tl] shows the density (black solid line) and tem- 
perature structure (black asterisks) of the g func- 
tion of vertical column density at radii R — 10, 20 and 
30 AU for the X-ray-|-EUV model with no pre-screening 
and Lx = 2 • 10^°erg/sec (FS0H2Lxl). The red dashed 
and the red dotted lines show the dust temperatures 
and densities, respectively, calculated by d'Alessio et al. 
(1998). 

The red asterisks indicate temperatures that are above 
the local escape temperature, and therefore the base of 
the flow is located by the red point at the highest col- 
umn density. Since a small error on the location of the 
flow could potentially translate into a large error on the 
density (and therefore photoevaporation rates) great care 
was taken to ensure that the Monte Carlo error is small 
in the region where the base of the flow is located. 

The gas temperatures show a small step at columns of 
~ 5 • lO'^^cm"^. This is not realistic and an artifact of 
our models which do not include a self-consistent treat- 
ment of molecular chemistry. We impose the edge of 
the "molecular zone" from the models of Glassgold et al. 
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Fig. 10. — Dependence of the total mass 
luminosity. 
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(2004), and assume that in the molecular zone fine struc- 
ture lines of neutrals are no longer efficient at cooling the 
gas (since all gas is molecular). This causes the tempera- 
ture to rise slightly as the gas and dust become thermally 
coupled. We stress that our temperature calculations in 
the deeper layers are not reliable and we postpone fur- 
ther discussion to future work when a chemical network 
will be introduced in the models. 

In the X-ray/EUV heated region the gas temperatures 
increase steeply and the gas densities respond with an 
initial steep decline, however both gas temperatures and 
densities flatten out at low column densities. This tem- 
perature behaviour is expected for a photoionised gas 
whose thermal balance in this region is dominated by 
heating by photoionisation and cooling by coUisionally 
excited line emission. The density behaviour is due to 
the fact that at high z's, where gravity is weak, the ex- 
ponential term in Equation 1 tends to unity and therefore 
the density at given height is then only proportional to 
the ratio of the gas temperature at the lower boundary 
point to the local gas temperature, i.e. the hydrostatic 
structure is roughly isobaric. 

This flattening of the densities at large heights demon- 
strates somewhat a failure of the HSE assumption. In 
reality, material in the upper disk atmosphere will be 
part of a flow and not bound to the disk itself. This 
should not have a significant effect on the photoevapo- 
ration rates, given that the isobaric layers are not in the 
line of sight between the source and the flow surface. 

5. COMPARISON WITH PREVIOUS WORK 

Previous studies of the thermochemical properties and 
photoevaporation of X-ray irradiated gas in protoplane- 
tary disks include ACP04, Glassgold et al. (2004, 2007), 
Ercolano et al. (2008, Paper I) and Gorti & HoUenbach 
(2009, GH09). 

ACP04 used cloudy (Ferland, 1996) (a one- 
dimensional photoionisation code) to construct a sim- 
ple 1+lD model to investigate the importance of X-ray 
photoevaporation for the dispersal of gaseous protoplan- 
etary disks around solar mass stars. Using the simple 
pCs approximation they obtained a surface-normalized 
mass loss rate E = 2.6 • 10~^^ g s~^cm~'^ at a radial 
distance of 19.1 AU. Although their model was rather 
idealised and they used slightly different input parame- 
ters for the central star and illuminating spectral shape 
and luminosity, the values of S they estimated is in 



the range of those obtained in this work (see Figure |4]), 
and corresponds to a total mass loss rate of roughly 
10~^ MQ/yr. ACP04 compared their estimate of E with 
that previously obtained for UV photoevaporation by 
HoUenbach et al. (1994), and found comparable val- 



ues (E 



UV 



3.88 • IQ-^^gs-^ 



at 6.7 AU), which 



lead them to the conclusion that X-ray photoevapora- 
tion does not play an important role. However, at a much 
larger radial distance of 19 AU — approximately the re- 
gion that Ercolano et al. (2008) found to dominate X-ray 
photoevaporation — ACP04 obtained similar values of S. 
Due to the R'^ dependence of the total mass loss rates, 
this translates to a photoevaporation rate nearly an or- 
der of magnitude higher than that at 6.7 AU. Therefore, 
while we agree with the numerical results of ACP04, we 
interpret the result differently and stress that our calcu- 
lations (and those of ACP04) both suggest that X-ray 
photoevaporation does play a significant role in disk dis- 
persal. 

Glassgold et al. (2007) calculated the thermochemical 
structure of a disk around a 0.5 Mq pre-main sequence 
star, irradiated by X-rays from a central source. The 
authors focused on the characterisation of the physical 
properties of the disk and the prediction of gas-phase 
diagnostics and do not give an estimate for X-ray photo- 
evaporation from their models. However from the model 
grids which were kindly made available to us by these au- 
thors, we flnd E = 2.2 ■IQ-^'^g 3-^ cm''^ at 20 AU. While 
this is higher than the estimates of this work and ACP04, 
it agrees very well with the values reported in Paper I. We 
note that Glassgold et al. (2007) and Paper I both used 
a fixed density structure for their calculations, and the 
higher photoevaporation rates are due to the fact that 
no obstruction by the puffed-up inner disk is accounted 
for in these models. 

Contrary to this work and to the results mentioned 
above, recent work by Gorti & HoUenbach (2009, GH09) 
predicts much smaller photoevaporation rates due to di- 
rect X-ray heating and concludes that X-rays are only 
of indirect importance for the dispersal of protoplane- 
tary disks around solar mass stars. Instead, their mod- 
els predict that X-ray ionization generally acts to "am- 
plify" FUV photoevaporation rates. While the authors 
state that their results are in agreement with those of 
ACP04, inspection of their Figure 5 (curve labeled X) 
shows that for their "X-ray only" model they estimate E 
of the order of 10~^^(7 s~^cto~^ at 20 AU, which is two 
orders of magnitude smaller than our and ACP04's esti- 
mates. We also note that ACP04 obtained their estimate 
of E 

than 10 AU as stated in the GH09 paper 

Comparison of complicated theoretical calculations can 
be difficult and the reason for disagreement may be 
buried in the codes themselves or in the physical assump- 
tions adopted. However we have identified at least one 
reason for the discrepancy between our calculations and 
those of GH09 which is the hardness of the illuminating 
X-ray spectrum. 

Indeed, we find that the soft-X-ray region is the most 
efficient at heating the gas in the photoevaporation layer. 
GH09 use a rather hard spectrum, which peaks at 2keV 
and is approximated by a power-law with Ahx /dE ^ E 
for O.K < 2 keV and dL^/dE - E^^ for 2 < < 10 



10 "^^gs ^cm ^ at a radius of 19.1 AU, rather 
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Fig. 11. — Gas temperature (black asterisks) and density struc- 
ture (solid line) for the unscreened X-ray+EUV model irradiated 
by a source at 2 R* (FS0H2Lxl). The red dashed line represents 
the dust temperature distribution of the d'Alessio (2001) model. 
The red asterisks mark the locations where the gas temperature 
exceeds the local escape temperature. 

keV. This is a significantly harder spectrum than that 
used by ACP04 (which peaked at 0.7keV), ourselves in 
Paper I (see Figure [T]) and by Glassgold et al. (2007, 
which peaked at IkeV). In this work we found that mod- 
els irradiated by a hard spectrum such as that produced 
by high column screens (e.g. FS22 models) did not yield 
significant photoevaporation. 

We have run a model using the GH09 spectrum and 
we obtained E ~ 6 • lO'^g s-^cm-^ at 20 AU (us- 
ing the pCs approximation) and a total mass loss rate 
of ~10~^° M^lyr. While a step in the right direction, 
the hardness of the illuminating spectrum alone docs not 
seem to be sufficient to explain the discrepancy and other 
factors must also come into play. 

6. LIMITATIONS OF CURRENT MODELS 

6.1. The estimation of photoevaporation rates from the 

models 

By far the biggest uncertainties concern our attempt 
to use these hydrostatic radiative equilibrium models to 
estimate the resulting mass loss rates. This now be- 
comes a critical problem given that our preliminary esti- 
mates suggest that, contrary to previous claims (ACP04, 
GH09), X-ray photoevaporation may be the major agent 
of disc dispersal. 

The uncertainties stem from the fact that we have de- 
rived the modified disc structure under the assumption 
of HSE, and although this probably gives a reasonable 
description of the density structure at small z, it is obvi- 
ously the case that a static description does not produce 
a corresponding mass flow rate. Instead, we have em- 



ployed the simple approach of deriving the hydrostatic 
structure and then, at the height where the disc mate- 
rial exceeds the local escape temperature, assume free 
sonic outflow from that point. A similar approach in 
the EUV case is less questionable for two reasons: i) the 
steep temperature discontinuity at the ionisation front 
and the near isothermal conditions above it means that 
the location of the launch surface is independent of de- 
tailed modeling and just depends on the regions of the 
disc where photoionised gas is unbound and ii) the den- 
sity at the launch surface depends only on conditions 
above the ionisation front and not on the structure of 
the underlying disc. Thus the fact that the disc below 
the launch surface is not strictly static has no effect on 
the resulting outflow rates. 

In the X-ray case, the temperature rises much more 
gradually with height: for example, for our model at 20 
AU, the temperature rises from roughly 200 to 5000 K 
over a region about 6.5 AU in height, over which the 
density of the hydrostatic structure falls by one order 
of magnitude. In reality, pressure gradients will cause 
the flow to be initiated below our nominal launch point, 
i.e. at a height where its temperature is less than the 
escape speed, and the density structure will be progres- 
sively modifled over the region where the flow is acceler- 
ated. It therefore becomes little more than an exercise 
in educated guesswork to assign a point for estimating 
the escaping mass flux - the lower down in the struc- 
ture one tries to do this, the more reliably the density 
structure will be known, but the less easy it is to assign 
a flow speed. At greater heights, it's easier to assign 
a flow speed (unbound gas expands at about the sound 
speed), but the estimate of local density becomes unreli- 
able, as the underlying density structure is now consid- 
erably modifled by the flow. 

It is tempting (see e.g. GH09) to propose a simple, 
physically motivated prescription that maps the struc- 
ture of the photo-heated region onto an analytic flow 
solution at large radius (for example, the isothermal 
Parker wind solution) . The problem with this is that al- 
though at large radii the flow will indeed approximate a 
quasi-spherical radial outflow, the transition between the 
nearly hydrostatic structure at the base and the Parker- 
like structure at large radius is quite different from a 
Parker wind - in particular, when the flow is non-radial, 
the interplay between terms in the momentum equation 
involving the divergence of the velocity, the gravitational 
acceleration and the centrifugal acceleration is quite dif- 
ferent from the non-rotating Parker case. The problem 
boils down to the fact that, until one knows the topol- 
ogy of the streamlines, one cannot compute the varia- 
tion of quantities along the streamlines - but that the 
steady state topology of the streamlines itself depends 
on the variation of quantities along the streamlines: in 
other words, it is a truly two dimensional problem. (See 
Begelman et al., 1983 for an attempt to parameterise 
the problem in terms of an ad hoc variation of the diver- 
gence of the velocity along streamlines). Our initial hy- 
drodynamic investigations of this problem (Alexander & 
Clarke in prep.) demonstrate that the streamline topol- 
ogy and the mapping between base density and mass flux 
is very dependent on how the base density varies with ra- 
dius in the disc. We therefore think that it is presently 
unwarranted to try and go beyond the simple estimates 
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employed here and that further progress in estimating 
photoevaporation rates will require a combined hydro- 
dynamics/radiative transfer approach. 

A further shortcoming is that a molecular chemistry 
network is still not included in our calculations. This 
poses a limitation on the accuracy of our gas temperature 
calculations in the lower layers where cooling by H2 and 
CO rovibrational and rotational lines may become impor- 
tant. For example at a radial distance of 1 AU and col- 
umn densities of 10^^-10^^ cm^^ Glassgold et al. (2004) 
find that CO rovibrational lines dominate disk cooling. 
CO rotational lines were also found to contribute some- 
what to the cooling at greater depths, although in these 
regions the thermal balance is dominated by dust-gas 
collisions (GN104). We are currently in the process of 
developing a chemical network to be included in future 
studies and will therefore be able to derive a more accu- 
rate thermal/hydrostatic structure at low temperatures 
{'^ 10^). It is questionable, however, how much such im- 
provements will affect our estimates of photoevaporation 
rates, given that gas in this temperature range is strongly 
bound at radii less than 100 AU. 

7. SUMMARY 

Our calculations of the structure of X-ray irradiated 
discs in low mass pre-main sequence stars imply total 
photoevaporation rates that, taken at face value, exceed 
the photoevaporation rates produced by extreme ultra- 
violet (EUV) radiation by an order of magnitude. This 
difference is attributable to the fact that a much more 
extended region of the disc is subject to X-ray photoe- 
vaporation than in the EUV case (i.e. significant evapo- 
ration out to ~ 40 AU, in contrast to the EUV case where 
little mass is lost beyond ~ 10 AU). The reason for this 
different photoevaporation profile is that the X-rays are 
able to penetrate columns of up to ^ lO^^cm"^ through 
the disc's puffed up, X-ray irradiated atmosphere and 
thus still produce significant mass flows when they im- 
pact the disc at several tens of AU. In the EUV case, by 
contrast, the regions of the disc from which the flow can 
be launched are shielded from direct irradiation by stel- 
lar photons and are instead subject only to recombina- 



Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2004, MNRAS, 
354, 71, ACP04 

Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2005, MNRAS, 
358, 283 

Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 
369, 216 

Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 
369 229 

Alexander, R. D. 2008, MNRAS, 391, L64 

Asplund, M., Grevesse, N., & Sauval, A. J. 2005, Cosmic 
Abundances as Records of Stellar Evolution and 
Nucleosynthesis, 336, 25 
Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 
Begelman, M. C, McKee, C. F., & Shields, G. A. 1983, ApJ, 271, 
70 

Bohren, C. F., & Huffman, D. R. 1983, New York: Wiley 1983, 
Clarke, C. J., Gendrin, A., & Sotomayor, M. 2001, MNRAS, 328, 
485 

D'Alessio, P., Canto, J., Calvet, N., & Lizano, S. 1998, ApJ, 500, 
411 

D'Alessio, P., Calvet, N., & Hartmann, L. 2001, ApJ, 553, 321 
D'Alessio, P. 2003, Revista Mexicana de Astronomia y Astrofisica 

Conference Series, 18, 14 
Drake, J.J., Ercolano, B., Flaccomio, E., Micela, G. 2009, ApJL, 

in press 

Ercolano, B., Barlow, M. J., Storey P. J., & Liu, X.-W. 2003a, 
MNRAS, 340, 1136 



tion photons from the bound region of the disc's ionised 
atmosphere. Since gas at lO^K is only bound out to 5 
AU, most of the EUV induced mass loss is restricted 
to regions modestly greater than this (HoUenbach et al., 
1994, Richhng & Yorke, 2000; Font et al, 2004). Our con- 
clusion that X-ray photoevaporation is more important 
than EUV photoevaporation is only strengthened when 
one considers two additional factors: i) the possibility 
of further attenuation of the EUV flux by absorption by 
any neutral flows (such as winds) close to the star and ii) 
the fact that we have unambiguous evidence that young 
stars remain strong X-ray sources even after disc dis- 
persal and thus are clearly available as a potential disc 
dispersal mechanism. 

However, we emphasise that our estimates of X-ray 
photoevaporation rates are much more uncertain than 
the EUV values, notwithstanding our careful treatment 
of two-dimensional radiative transfer and our derivation 
of the puffed up density structure of the irradiated disc 
(see discussion in Section 6). These uncertainties were 
of academic interest as long as it seemed likely that X- 
ray photoevaporation was in any case a minor agent for 
disc dispersal. We argue that future radiation hydrody- 
namical simulations are required in order to make firm 
predictions about the role of X-ray photoevaporation in 
disc dispersal. 
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